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We investigate the influence of the Jahn- Teller distortion and a direct antiferromagnetic moment 
coupling as extensions to a two-band Kondo lattice model for the magnetic and electronic properties 
of manganites. Those are calculated self-consistently via an interpolating self-energy model and a 
modifled RKKY technique using finite Hund coupling and quantum spins. We found that both ef- 
fects are essential to achieve realistic Curie temperatures if we regard intraband Coulomb repulsion. 
Using reliable model parameters we got results which are in very good agreement with experimental 
data in the whole ferromagnetic doping range. In the calculated phase diagram there are ferromag- 
netic metal to paramagnetic insulator transitions, accompanied by a Colossal Magnetoresistance 
(CMR) behavior. To improve the comparability of the measured behavior of the resistivity with the 
calculated one, we have to switch on interband Coulomb correlations. 



I. INTRODUCTION 



Manganites like Lai^ajCa^jMnOa or Lai_2:Sr2:Mn03 
obtained a lot of interest in the last decades, especially 
after the discovery of the resistivity's dependence on the 
magnetic field, the colossal magnetoresistance (CMR). 
The detection of huge changes of the resistance in the 90 's 
led to an acceleration of the research in this topic. Ex- 
perimentalists found very rich phase diagrams containing 
para- and ferromagnetic insulating and metallic regions 
as well as different kinds of antiferromagnetism, which 
can be accompanied by an ordering of charge or of or- 
bitals. Most of these phases can also be identified in 
theoretical calculations. But even after that long time 
of research there are still many unresolved issues. The 
origin of the CMR is determined as a competition of dif- 
ferent phases, but the exact connections are unknown. A 
large variety of phases is an evidence of the complexity of 
those systems. Thus it is important to create a solvable 
model, which contains as many as possible effects that 
can appear in manganites. 

The electronic and magnetic features are mainly due to 
the 3(i-electrons of the Mn'^^ and the Mn'^^ ions, whose 
ratio is defined by the doping rate x. Both kinds of ions 
have fully occupied spin-up t2g-\eve\s which provide a lo- 
calized spin of S* = I , but only the Mn'^+ ions also have 
an itinerant e^-electron. The main features of this system 
should be covered by the Kondo lattice model (KLM) or 
double exchange model. But already ten years ago it 
became clear that this model needs to be extended to 
describe manganites in a proper way. First of all Millis 
et al. argued that the Jahn- Teller effect (JTE), which 
can split up the e^-orbitals, should be important at low 
and intermediate doping rangel. Likewise it is stated 
in some recent publications that a direct antiferromag- 
netic coupling between the localized spins, caused by a 
superexchange via the oxygen orbitals, is necessary to 
describe the competition between the different magnetic 
phases . Besides the necessity for the different phases, 
it also seems to be important to increase the CMR in 



theoretical calculations^i^i^i^ to be in better accordance 
with the measured data. 



II. MODEL 

In this article we want to describe the manganites in 
the whole ferromagnetic regime. Therefore we use a 
two-band KLM with finite Hund coupling. We extend 
this model by terms, which describe the JTE, the su- 
perexchange and the Coulomb interaction. In that gen- 
eral model most of the physics of the manganites should 
be incorporated. The phonons of the JT modes are 
treated classically, but the spins quantum mechanically. 
The model is solved approximately, but self-consistently. 
These extensions will appear essential to get the cor- 
rect doping dependence of the Curie temperature and to 
achieve metal-insulator transitions simultaneously with 
the breakdown of the ferromagnetic order. Therefore we 
study the Hamiltonian: 

n = Hs + H,a + Hu + HAF + HjT (1) 

+ ^ V^^, n^aariia'a' + J AF Sj ■ S., ~ 

i,","! {i,j} 

i 

^ial annihilator (creator) of an electron with 

spin a in the orbital a corresponding to the d^2_j,2 or 

d^z^-r^ orbital at site i and riiaa = cl^^Ciaa- Jh is 
the Hund coupling between the itinerant spin cria and 
the localized S.^. The Hubbard parameter can represent 
the fact that there is only intraband Coulomb repulsion 
(i.e. C/"^? — 0) or that there also is an extra in- 
terband interaction {V^^ 7^ 0). Jaf corresponds to 
the direct superexchange. Q2,d,i are special JT modes 
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which couple with the pseudospin operators T^^'^, e.g. 

= ^J2cri'^i'ya=+i - niaa=-i)- A similar model is 
used e.g. in Ref.->i. Possible other effects, like a hy- 
bridization between the e^-bands, will not be considered 
explicitly. We will have a short discussion about that in 
the appendix El 



III. APPROXIMATION METHODS 

A. Electronic part 

To solve the resulting equations we need some sim- 
plifications. The JT part is treated in a mean-field 
approximation^. Assuming translational invariance, in- 
troducing spherical coordinates Qs = Qsin0, Q2 = 
Qcoad and using "dressed" operators c[^l, which are 
now corresponding to the higher and lower JT orbital, 
yields 



n = 0.1 



n = 0.5 



0.96 



kacr 

Here is Za^±i = =bl and 



(An) =-Y,z^{ 



(2) 



(3) 



As one can see (An) is dependent on the occupation 
numbers {na„) and therefore has to be calculated self- 
consistently. The dressed operators are 

— %Q 11 1 ^ ^ 

Ciaa=~l = e ' (cos -Ci£r,3z2-r2 + sin -Cio.^a,2_j^2) (4) 

Ci<TQ=+i = e*''/^(sin^Cic,^3^2_^2 -cos^Cio.^2.2_j,2) 

and are now a superposition of the original e^-states. We 

simplify the notation by setting c^^^ c^^ct- "^^^ formu- 
lation of HjT in is possible because of the restriction 
to non-cooperative effects. Of course, the JT-distortions 
are cooperative in reality, which is e.g. important for or- 
bital ordering. But we restrict on the non-cooperative 
distortions, so that we can solve the resulting equations. 
This suitable form of Hjt allows us to combine it easily 
with the kinetic part 



kacr 



(5) 



containing the dispersion e(k) of the free band for simple 
cubic structure. In this form of Hjt one can easily see 
that the JTE can split up the bands for (An) > (band 
Jahn- Teller effect). That is why we will now refer to both 
bands as JT bands, which are denoted by the index a. 
If we first leave out the antiferromagnetic coupling, we 
remain with a correlated KLM, 
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Figure 1: (Color online) Differences of the QDOS at ferromag- 
netic saturation (T = GK) if one only has intraband Coulomb 
interactions or additional interband repulsions. The main dis- 
crepancies are at higher occupation number n. Especially a 
complete filling of the bands is only possible with extra inter- 
band repulsion (Mott-Hubbard- insulator). The upper Hub- 
bard bands are shifted to infinity. J = leV, W — 3eV, g = 
0, Jaf = 



The related Green's function can be approximated within 
an interpolating self-energy approach (ISA)^. In this ap- 
proach we get a self-energy, which fulfills practically all 
limiting cases of the KLM. The exactly solvable cases 
of the ferromagnetically saturated semiconductor (mag- 
netic polaron) and the atomic limit are reproducible by 
the ISA. Also a second order perturbation theory is in- 
cluded. These cases are interpolated by the use of rig- 
orous high-energy expansion. If we postulate that this 
approximation is valid also between the limiting cases, 
we get a theory that should be reliable for all tempera- 
tures, band occupations and couplings Jh- 

Within the ISA, the Hubbard term was handled in an 
effective medium approach, so it influences the spectral 
weight and the width of the bands. But of course it will 
also give an additional splitting of the bands into differ- 
ent Hubbard sub-bands. We will call all bands (distin- 
guished by the indices {a, a) coming from the JTE and 
the spin) which are not affected by the Coulomb repulsion 
lower Hubbard bands and the other ones upper Hubbard 
(sub-)bands. The lower and the upper sub-bands have 
an energy difference of the magnitude of the Hubbard 
parameter U. For manganites the value of the Coulomb 
repulsion is much greater than all the other parameters 
so we can choose in a good approximation U —^ 00. This 
means the upper Hubbard bands will be shifted to infi- 
nite energy and will therefore never be occupied. Thus 
we do not have to consider these bands explicitly. We 
get the one particle Green's function 



,Ck, 



vISA 



= h 



la 



i? + M-T„,(k)-I]^SAj^) 



It is 



1 - (no 

intraband 



(7) 



(8) 



interband 



n' = H 



JT 



H, 



(6) 



the spectral weight of each band coming from the Hub- 
bard term, S{^^(i?) the self-energy in the ISA formalism 
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and TaaO^) = Zag'^i^n) +7ccre(k) the centers of gravity. 

Here the represent the probabihty that there is 
no repulsion partner for an electron on the same site. 
Thus the sum of the 70,0- , which is equivalent to the max- 
imum occupation of the lower Hubbard bands, for the 
mfraband repulsion is greater than for mierband inter- 
action. Figure [1] shows the influence of the -faa on the 
quasi-particle density of state (QDOS). As can be seen 
a complete filling of the lower Hubbard bands for n — > 1 
is only possible with the adding of interband repulsion. 
This complete filling correlates with the generation of a 
Mott-Hubbard-insulator. Thus this property of the man- 
ganites for n = 1 is automatically fulfilled within the ISA 
including an extra interband repulsion. 

The explicit structure of the self-energy is 



viSA 



(9) 



H' 



1 — hJHGf^_^{E — ^ZaJnXa^-a) 



where 



1 ^ao" 



(10) 



and = S'f +iz„Sf . With the spectral theorem we can 
calculate important terms like the correlation functions 



nNJ 



f dEf^{E) 



x[S-r,„(k)]ImGw(i;-^) (11) 
the mean occupation values 

M = [ dE f^iE)lmGk^^{E - ^i) (12) 

and the occupation difference 

(An) = ^ {{na= + l.a) - {na=-l.cr)) 

(7 

self-consistently. The only remaining model parameters 
are the spin S, Jh, the JT coupling g, the bandwidth W ^ 
the total density n = \ — x (corresponding to /i) and the 
magnetization (S*^). 

There are some important changes in the QDOS of the 
ISA compared to e.g. mean-field calculations. The often 
used neglection of minority spins by setting Jh ^ 00 con- 
tradicts for example the KLM's exactly solvable atomic 



limit. Even for large Jh there is a spectral weight of the 
spin-down peaks of the QDOS, so that a large Hund's 
coupling alone cannot prevent double occupation. Thus 
we get a finite QDOS of the spin-down electrons at lower 
energies (Fig. [T]). This originates from scattering pro- 
cesses of the spin-down electrons. Those can do a spin- 
fiip while emitting a magnon and so become a spin-up 
electron. This can only be done at an energy interval 
with a finite number of spin-up states, due to the low 
magnon energies. That means we have a finite occupa- 
tion number of minority spin electrons even for a large 
Hund's coupling. Furthermore there is a peak in the spin- 
up spectrum at higher energies which is connected to a 
magnetic polaron. Details can be found in the original 
paper—. 



B. Magnetic part 

To calculate (S^) in a self-consistent way, too, we use 
another technique. In the modified RKKY formalism 
(mRKKY) we try to map the Hamiltonian of the KLM 
onto an effective Heisenberg mode l^°i^^ . That is done 
by averaging out the electronic degrees of freedom. This 
yields 



cka' ; (13) 



-icrcr'kq 



v(«) 



means averaging while treating the spins as 
numbers. We can now construct the equation of mo- 
tion (EOM) for the Green's function G'£fk+q(-£^) = 



This yields 



G-^^+^{E)^S^oS..'G'^\e) 



i(k-k')Ri^(0) 



2hN 



E 



(14) 



Gk {E)(Si ■ <t) o-o-"Gk',k+q 



+ e 



-.(k'-(k+c))R.(5.^V(^)(S^ • tU'-g'^Ue) 



As a simplest approximation we can replace the full 
Green's functions r.h.s. by free Green's functions, 

G^i^k+qiE) Gj^+q(i^)'5k',k+q'5o-',o-" 

G^uM'iE) - G("^(i?)<5k',k<5.,.< (15) 
Then we use the spectral theorem and 

' '^2 W -Si 

to get effective exchange parameters 

j.y.O. l^,.,.y^ /-(Kk + q))-/_(.(k)) 
^(q) - 2^///^ 6(k-Hq)-e(k) 
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Figure 2: (Color online) Tc- Jfr-dependence in the mRKKY 
formalism, left: For higher densities n a Jh > Jh is needed, 
to get a finite Tc. The regions where the mRKKY corresponds 
to conventional RKKY and the double exchange formalism 
are marked. W = leV right: With increasing bandwidth W 
a higher Tc is achieved, but a larger Jh is needed to get a 
finite Curie temperature, n = 0.95 both: g = 0, Jaf ~ 0, 
intraband repulsion 
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Figure 3: (Color online) Influence of the JT-coupling on Tc. 
Beyond a critical gc (marked by arrows) the JT-splitting of 
the fig-bands starts and then Tc is changed. If there is only 
intraband Coulomb repulsion Tc is everytime decreased, but 
can be increased if we switch on an additional interband re- 
pulsion. Jh = 2eV, W = leV 



for an effective Heisenberg model 



4 - ^E-^(q)^""^''-''^ 



This is the result of the conventional RKKY (cRKKY) 
which is originally derived by perturbation theory for low 
Jh- But we can go one step further and replace the upper 
Green's functions by the Green's function we derived in 



where J^^ has only next-neighbor elements with the 
magnitude Jaf- It is important to recognize the dif- 

(17) 

ferences between both couplings. comes from the 
electronic properties and is a long range interaction (w 
5th-30th neighbor) and J^^ is constant and short range 
(next neighbor). Thus the adding of J.^^^ is more than 
just a renormahzation of the energy scale. We can now 
directly calculate the magnetization (S^) with the solu- 
tion by Callerti^: 



^k+q iE)5ii' M+qSa' ,cr" 



(18) 



After the solution of the EOM, we get effective exchange 
integrals 



i(q) 



+ 00 -1 

dEf_{E) 



kacT 



{E) (19) 



which define the effective Heisenberg term. These ex- 
change integrals now contain the electronic correlations 
due to 



^k,k+q 



{E) - G 



(0) 



ka(T(-^)G'k+qQCT(^) + 

h 



(20) 



E + ^i-T^,{y) 



where G^^^{E) is the Green's function ([7]). We have 
now included the band index a which we first left out for 
simplicity. The mRKKY covers the cRKKY for low Jh 
{Tc ^ J]f) and the double exchange for large couplings 
{Tc ~ const) (Fig. [2]). After the mapping onto an effective 
Heisenberg model, it is easy to add the antiferromagnetic 
part by just summing both exchange integrals 



■I 



AF 



(21) 



{1 + S + + {S ~ ip){l + ip) 



2S+1 



-E 



(1 + <^)2S+1 
1 



S(q) = 2^(5^) ( Jo - J(q)) - gj/iis^o 



q 



(22) 



Furthermore one can derive equation (jBip for the direct 
calculation of the Curie temperature as it is shown in the 
appendix. 

A typical feature of the mRKKY for larger n is the 
appearance of a critical Jj., which is needed to get a finite 
Tc- This Jc can become very large, especially for a larger 
value of the free bandwidth W. It can even achieve the 
magnitude of the large Jh of the manganites (Fig. [2]) . 

The whole model is now self-consistent with respect to 
the electronic and magnetic properties. 



IV. RESULTS 

We will focus on the model's electronic and magnetic 
properties, like the resistivity, the JT-splitting and the 
Curie temperature. As we will see there is a great dif- 
ference between having only intraband repulsion or hav- 
ing additional interband correlations. Interesting, for ex- 
ample, is the influence of the JT-splitting on the Curie 
temperature (Fig. [3]). In our self-consistent calculation 
we need a critical which was also found in other 
works^'liSi^ to split up the bands . This splitting will 
usually decrease Tc, but with extra interband repulsion 
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and larger occupation number n (n > 0.8, which means 
half-filling in this case) T^, can also be increased. Thus we 
will observe the respective interaction in different parts. 

Experimental results show that the ferromagnetic 
regime is in the doping rate 0.1 ^ a; < 0.5 with a maxi- 
mum of Tc at a; ~ q 322,23 ^pj^g observed behavior of the 
Curie temperature is not reproducible in the pure two- 
band KLM, because it causes a finite Tc for high and 
low doping rates. Therefore it is necessary to include 
other effects. In our case these are the JTE and the su- 
perexchange, incorporated by a direct antiferromagnetic 
exchange Jaf ■ Furthermore at higher electron densities a 
CMR is observed accompanied by a FM to PI transition. 
To achieve a large CMR and to reproduce the insulating 
behavior for n — > 1 we will need interband Coulomb re- 
pulsion. 

Several suggestions were made to explain the CMR. Most 
are based on a competition between different phases (e.g. 
Ref.^). In Ref.^ this is described as a result of the com- 
peting tendencies to form a charge localized or a ferro- 
magnetic phase. This leads to a peak of the resistivity at 
Tc- A further theory is based on the current-carrier den- 
sity collapse induced by the emerging of bipolaronsi^J^. 
Those bipolarons can be broken up below Tc for a ex- 
change JhS which is large enough. We will see in our 
model that this peak of the resistivity at Tc can be the 
consequence of a lattice distortion (intraband Coulomb 
repulsion) or of a drastic change of the spectral weight of 
the quasi-particle density of states (interband Coulomb 
repulsion) . 

It is not in the scope of the paper to investigate an- 
tiferromagnetic phases or orbital ordering. Thus we will 
only focus on ferro- or paramagnetic phases. 



A. Intraband Coulomb repulsion 

First we will only use intraband Coulomb repulsion. 
That means we use ja^y = 1 — (na.-a) in 0. Figure [3] 
shows how the JTE and the antiferromagnetic coupling 
act on the Curie temperature in this case. An increase 
of the according couplings reduces Tc- The impact on 
Tc of the JTE appears at low and of the AF coupling at 
high doping rates. This different scope of each couphng 
is physically consistent, because at high electron density 
(that means small x) there are more active Mn'^"'" ions, 
so the JTE should be supported, which is shown by mea- 
surements, too2i. A reduction of Tc by the JTE has also 
be found by other authors, e.g. Refill. The drop of 
the Curie temperature at small x can also be reproduced 
without the JTE? and/or if we include a large interband 
Hubbard interaction in our model. 

Jaf was found to stabilize the antiferromagnetic CE- 
phase for x > 0.5^^, so it seems plausible that the influ- 
ence of this coupling on Tc should be noticed stronger in 
this region. The importance of Jaf in connection with 
the RKKY- mechanism is even more obvious. No matter 
whether we use the cRKKY or the mRKKY formalism. 



AF coupling JT coupling 




Figure 4: (Color online) Effect of the different couplings on 
the Curie temperatures at intraband Coulomb repulsion, left: 
Jaf acts on Tc mainly at higher doping rates and can re- 
duce Tc to zero in this regime. This is not possible in the 
pure KLM. right: The JT-coupling acts at lower doping rates 
and can change Tc only in a special range r™'"(2;)c < Tc < 
T^^^{x), in contrast to Jaf- The Curie temperature will not 
be reduced for g > Qoo any more (comp. Fig. [S}. The up- 
per Tc-curve on each side matches the normal two-band KLM 
without AF and JT coupling, left: g = rigtit: Jaf ~ both: 
Jh = 2eV, W = leV 



there is always a finite Curie temperature for small elec- 
tron densities, i.e. high x, for Jh > (compare Refj^i 
and Fig. I2l4p . Thus we need a non-RKKY-like effect to 
reduce Tc in this region. It should favor antiferromag- 
netism, which is observed for higher doping rates. To 
the knowledge of the authors, there is no other plausible 
effect, which reduces Tc at higher x and shifts its maxi- 
mum to lower doping rates and which can be added to the 
two-band KLM. Calculations without such an extension 
can achieve good results for Tc at lower doping rates, but 
they will miss the characteristical decrease above x « 0.3. 

We can now try to calculate realistic values of Tc. 
Therefore we choose = | and fix the Hund cou- 
pling to Jh = 2eV, which is a typical value and con- 
firmed (with large error bars) by photoemission measure- 
ments. For the free bandwidth we choose W = 2eV 
for Lai_j;Ca2,Mn03 and Ndi_2:Sr^Mn03, which is esti- 
mated from ab-initio calculations for Lai-^CaajMnOs — . 
Since Lai_xSra;Mn03 is regarded to have a larger band- 
width, we use W = 3eV. That corresponds to a hop- 
ping t — 0.16(0.25)ey for the simple cubic model den- 
sity. The order of magnitude of the JT and AF cou- 
pling can only be determined roughly. In Ref.* g is 
guessed g = 1 . . . l.G^/kjrt, which means in our model 
with kjT — 1 that gt=o.i6eV = 0.4 . . . 0.64\/eT/. Jaf is 
considered to be in the range of Jaf — 0.01 . . . O.li^i^l => 
Jaf ~ 2 . . . 2QmeV. These parameters will be fitted to 
the experiment for each material. That means we choose 
a set of these parameters and calculate full magnetization 
curves {Sz){T) for each x and get Tc from that curves. 
It is necessary to calculate the full curves because of the 
temperature dependence of the strain and the mutual 
influence of the JTE and the magnetization (discussed 
later). If one calculates Tc directly by using (jBip one has 
to set (Sz) = and one has no influence of the magneti- 
zation on the JT splitting any more. 
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Figure 5: (Color online) The calculated Curie temperatures 
in comparison with the experimental data for the correspond- 
ing materials in the ferromagnetic regime^>2^. Lines are a 
guidance for the eyes. The doping range where a JT split- 
ting occurs is marked. For Lai-xSr^^MnOs there are typical 
curves without the JTE or Jaf but fitting the other param- 
eters. One can see that both effects are needed to get the 
typical shape of the Tc-curves. Parameters: Lai-^^Ca^MnOa 
W = 2eV, Jaf = 5AmeV,g = 0A28VW Ndi_,Sr,Mn03 
W = 2eV,JAF = 5.0meV,g = 0.435^eF Lai_^Sr^Mn03 
W = 3eV, Jaf = 3.7meV, g = 0.436\/iF all: Jh = 2eV 
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Figure 6: (Color online) left: Comparison of the calculated 
values for the strain and (Sz) for Lai_a;Ca:rMn03 at x = 0.15 
and the measured strain^*. The magnetization first lowers 
the JT distortion and then breaks down at the rising of the 
strain at Tc. right: Phase diagram for the parameters of 
Lai-ajSr^iMnOs in Fig. \5\ M metal, / insulator, F(P) ferro- 
(para-) magnet. At the M-I transition occurs a CMR effect of 
a maximum of 170% at a; = 0.1. 



As Fig. [5] shows it is possible to have an excellent agree- 
ment with the experimental data within our theory. The 
varied parameters g and Jaf stay in the estimated range 
and play a crucial role to achieve the conformity of the 
curves of the Curie temperatures. Firstly, which is most 
important for an approximative theory, our method can 
reproduce the right trends according to Tc and secondly 
it can also give qualitatively right values in contrast to 
e.g. mean-field treatments. Thus it could be applicable 
to be combined with ab initio calculations to derive self- 
consistently input parameters like Jh- 

In the low doping region exists a finite JT splitting of 
the two JT bands of the order Ejt = 2g^(An). Near 
the occupation difference (An) (equivalent to the strain 
Q — g{An)) increases, which is in qualitative agreement 
compared to^^ (Fig. [6]). Even though the transition tem- 
perature of the distortion is not far away from the ex- 
perimental result, the main focus should be on the qual- 
itative behavior, because of the simple treatment of the 



JT-term. 

The increase of the strain is accompanied by a drastic 
decrease of the magnetization. That is the cause why we 
could not use formula (|B1[) to calculate Tc- If we are in 
the critical regime of the parameter g the JT splitting and 
the magnetization have a large impact on each other'^'^^. 
At intraband repulsion both effects suppress each other. 
That means the larger the magnetization (5*2) the lower 
the JT splitting Ejt and vice versa. When (S^) goes 
down because of the rising temperature it cannot have 
that much effect on the JT splitting any more. Therefore 
the splitting of the JT bands becomes stronger. With this 
increase of Ejt it will now suppress the magnetization 
even more and will lower it in addition to the lowering of 
the temperature effect. When this self-energizing effect 
exceeds a special value, we have a first-order transition 
of the magnetization. One can also see this behavior in 
the QDOS (Fig. [T]) where the different bands split more 
for T = Tc- The finite QDOS for the spin-down electrons 
at lower energies is a quantum mechanical effect due to 
scatterings. 

If we calculate the resistivity according to (|C2[) and 
define an insulator/metal via the temperature behavior, 
we can get the phase diagram of Fig. [6l There are metal- 
insulator transitions at lower doping rates, which are ac- 
companied by a change of the magnetic phase. In Fig. [5] 
(left) the explicit resistivity curves are plotted. One sees 
the jump of pci at Tc for low doping rates. The reason 
for that is the behavior of the JT-splitting (comp. Fig. 
El left). At Tc the bands split which means a lower total 
DOS at the Fermi energy (comp. Fig. [7]) and therefore a 
decrease of the conductivity. With rising temperature the 
JT-splitting becomes smaller and the bands get more and 
more overlap again. Thus the resistivity decreases with 
temperature resulting in an insulating behavior. With 
the breakdown of the JT splitting we get the normal 
metallic phase again. Though the simultaneous FM/PI 
transition leads to a CMR behavior, the resistivity jump 
is too small compared to measured results. 

Antiferromagnetic phases and phase separation are 
very important, but were not part of this work, because 
it would enormously expand the complexity of the prob- 
lem. Those were found in Refi^ for the KLM using the 
same self-energy. The inclusion of the AF-phases within 
our whole model will be left for later investigations. 



B. Interband Coulomb repulsion 

With intraband repulsion we achieved very good re- 
sults concerning the Curie temperatures and the JT split- 
ting. Even though there arc important phases and phase 
transitions (FM-PI), the ferromagnetic insulating (FI) 
phase is missing and, of course, the many antiferromag- 
netic and charge /orbital ordering phases are not there, 
either. While we will not consider antiferromagnetism 
and orbital/ charge ordering, we will try to get the FI 
phase in our model. Such a phase is found for low dop- 
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Figure 7: (Color online) The quasi-particle DOS and 
fermi function (blue solid line) for the parameters of 
Lai-^Caj^MnOa in Fig. [5] at 2: = 0.15, below and at Tc = 
175K. Those bands originate from the combined states of 
Cg-orbitals shown in (|4]). One sees that the splitting (An) 
between the lower (solid line) and the upper (dashed line) JT 
band rises at Tc and the same occurs for the strain Q — g{An) 
in accordance to the experiment^*. For a large Jh there is a 
finite occupation of the spin-down band for T < Tc, too. 



ing rates, e.g. in Ref.— . This phase cannot be repro- 
duced in the intraband model. These measurements also 
show that the resistivities decrease with increasing dop- 
ing rates. That is not reproduced in the intraband treat- 
ment, thus we will now add the interband Coulomb in- 
teraction. 

Actually we can achieve a ferromagnetic phase with 
insulating behavior of the resistivity (Fig. [8|). Addition- 
ally we have a better comparability of the theoretical and 
experimental results for the resistivity in the whole ferro- 
magnetic regime, too. In particular the system becomes 
an insulator at a; ^ 0, the Mott-insulator as described in 
Fig[TJ But the ferromagnetism breaks down in this limit 
(Tc = OK) which indicates that other phases are impor- 
tant for x — > 0. This is indeed the case, as it is widely 
known from the experiments showing antiferromagnetic 
or spin glass phases, which we are not treating in this 
work. That our approximative self-energy T}J'^{E) actu- 
ally leads to antiferromagnetism at a; ^ was shown by 
Hennig^S.. 

Furthermore we can get a high CMR effect (Fig. [9|) 
which was very much smaller while using only intraband 
repulsion. In this case the origin of the resistivity jump 
is different from that which occurred in the intraband 
calculations. The latter originated from the JTE and 
the increasing band splitting. Now, in the interband cal- 
culations, it comes from a drastic change of the DOS 
in the bands themselves. Figure [TO] shows that at it 
comes to a large rearranging of the DOS at the Fermi en- 
ergy if no external magnetic field is applied. This leads 
to a drastic reduction of the spin-up DOS which cannot 
be compensated by the increase of the spin-down DOS. 
Therefore the resistivity is jumping to a higher value. At 
higher temperatures above Tc there is an increase of the 
conductivity due to the softening of the Fermi function. 
The fast changing of the DOS can be slowed down in the 
presence of an external magnetic field and leads to a de- 
lay of the increase of pci- That means we have a CMR 
effect. 

Quantitatively, the resistivity should vary over more 
orders of magnitude, according to the experiment. We 
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Figure 8: (Color online) left: Calculated resistivities for the 
parameters of Lai-^^SrajMnOs in Fig. O intraband repul- 
sion right: Calculated resistivities with extra int erband re- 
pulsion: W = 2eV, J = SeV, Jaf = 0, 5 = O.SVeV The values 
of Tc are marked by solid arrows and the critical tempera- 
tures of the JTE by dashed ones. With extra interband re- 
pulsion the curves now show qualitative agreement with the 
experimen1i22 and especially a ferromagnetic insulating phase 
can be achieved at lower doping rates. 




100 
T[K] 



Figure 9: (color online) The resistivity for different exter- 
nal magnetic fields. With extra interband Hubbard-repulsion 
CMR-behavior with Ap = p(h=o)-p(h ) ^^^^ ^qq^^ 

' P(«max) 

found. X = 0.1, g = 0^/eV, W = SeV, J = 3eV 



also have to pay for the better agreement of the resistiv- 
ity with a worse behavior of Tc. The Curie temperature 
is much more suppressed if we use interband repulsion 
and we cannot get the typical shape of a parabola for the 
Tc-n-curves. Such a strong influence on Tc is probably 
due to the effective medium treatment of the Hubbard 
part. This treatment cannot contain all many-body cor- 
relation effects. Thus for the calculation of Tc we should 
include the Coulomb interaction in a more subtle way, 
which is no simple task, of course. 



V. SUMMARY 

We investigated manganite systems with a two-band 
KLM, which was extended by terms that represent the 
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Figure 10: (color online) The quasi-particle DOS (solid lines) 
for parameters of Fig. [9] at different temperatures. Only the 
lower part of the QDOS closed to the Fermi edge (~ —2.15eV , 
dashed line) is shown, left: {B = OT) The QDOS at the 
Fermi energy is changing very drastically at Tc- It develops a 
lack of spin-up electrons which cannot be compensated by the 
spin-down electrons and creates the jump of the resistivity in 
Fig. [5] After Tc the QDOS is not changing very much and 
pel lowers because of the softening of the Fermi edge, right: 
(B = ler) The QDOS is not changing very much at all. Thus 
it exists just a small change of pci- 




Figure 11: (color online) Resistivy curves for different hy- 
bridization values V at two doping rates. Major discrepan- 
cies occur only if V is much greater than the normal intra- 
band hopping t = O.lGeV. Parameters: W = 2eV, Jh = 
3eV,JAF = 0,g — 0.5\/eV, interband repulsion (comp. Fig. 



properly and none of the parts seems to be negligible. 

We are grateful to G.G. Reddy for fruitful discussions 
during his stay in Berlin. 



Coulomb correlations, the JTE and the superexchange. 
Within this model we calculated the electronic and mag- 
netic properties self-consistently by the use of an in- 
terpolating self-energy approach and a modified RKKY 
method. Because of the use of full single-particle Green's 
functions, this method contains more many-body inter- 
actions than the conventional RKKY. Therefore it gives 
reliable results even for larger Jh- 

With this formalism it was possible to calculate Curie 
temperatures, which are in very good agreement to ex- 
perimental measurements in the total ferromagnetic dop- 
ing range, if we use intraband Coulomb repulsion. We 
have shown that therefore the additional terms to the 
KLM are essential to achieve these results. We found a 
phase diagram with FM-PI transitions where a CMR ef- 
fect occurs and the JT distortion behaves qualitatively 
like in the measurements. But we have seen in our model 
that the neglection of interband Hubbard correlations 
will lead to an incorrect behavior of the resistivity. With 
the introduction of such interactions we get qualitatively 
correct results. A disadvantage of our effective medium 
treatment of the Hubbard part is, that it has a too strong 
influence on the magnetic properties, especially Tc. Pos- 
sibly, a treatment of the Hubbard term that has a better 
inclusion of many-body correlations could correct these 
discrepancies. 

Most of the research on manganites is done with the 
simplification Jh oo, which means effectively the ne- 
glection of the minority spins. But as can be seen in Fig. 
[71 there exists also for T < Tc & finite occupation of the 
spin-down band. This comes from scattering processes 
of the spin-down with the spin-up electrons accompanied 
by magnon emission or absorption, which is a result of 
the quantum mechanical treatment. Thus the infiuence 
of the spin-down electrons can not be neglected, just by 
the assumption of infinite Hund's coupling. 

In principle our model contains all ingredients to de- 
scribe the para-/ferromagnetic phase of the manganites 



Appendix A: HYBRIDIZATION 

Hopping between different bands seems to be impor- 
tant in manganites. To recognize which effect such a hy- 
bridization could have in our model, we use an effective 
medium approach: 

kcra 

• • ■ + E^C^-"Ck.-a (Al) 

kcTQ 

Here the first part is the original Hamiltonian ([1]) , repre- 
sented by the approximative ISA self-energy (jlOp contain- 
ing the correlation effects. Due to the adding of Ejf^(i?) 
the electrons can be treated like free electrons in an ef- 
fective medium. The second one describes the interband 
hopping with the hybridization V . After solving the ac- 
cording EOM we can investigate the possible changings 
due to this new term. 

There are no big changes of the resistivity curves for 
hybridization values at the order of magnitude of the in- 
traband hopping {t « Q.2eV), as can be seen in Fig. [TT] 
Only if V exceeds this range it can come to essential 
modifications. This can be also shown at the density of 
states (Fig. [T2|) . The formation of a hybridization gap 
occurs only for larger V . For V <t the QDOS is almost 
unchanged. 

Actually those bands we have investigated do not come 
from pure d^z'^-r'^ or dx2-y2 eg- orbitals, but from the 
states shown in (jH). That means the resulting bands do 
not belong to one of those e^-orbitals, but already con- 
tain some mixing. 
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Figure 12: (color online) Quasiparticle DOS for different hy- 
bridization values V. Lower JT-band is represented by the 
(black) solid line, upper by the (red) dashed line and the 
Fermi function by a (blue) solid one. Only the part at lower 
energies closed to the fermi edge is shown. For hybridizations 
in the near of the intraband hopping t = O.lGeV there are only 
small changes. Parameters like those in Fig. [Tl]for a; = 0.15 
and T = 50K, interband Coulomb repulsion. 



Appendix B: DIRECT CALCULATION OF Tc 



parameter constellations. Thus (|B1|) has to be handled 
carefully. 



Appendix C: RESISTIVITY 

To calculate the electrical conductivity tensor, we can 
use the Kubo forniula^^ and get a current-current corre- 
lation function 

(7«/3(S) = F f d\ f (j'5(0)j"(i + iA;i))e^(^+*^)* dt 



For the special case e(k) = ^(^k) and Wa(k) — 
■|9fe^e(k) — — k), which holds for the simple cu- 
bic structure, this can be simplified to a formula, which 
only contains one-particle Green's functionsi^iil. With 
the definition of a transport function, e.g. in x-dircction. 



By setting {Sz) 0+ in equation one can get an 
explicit formula for the Curie temperature 



k^T^^ -hS{S + 1) 



V 



-E 



1 



J(0) - J(q) 



(Bl) 

This formula (|B1[) is exact, but it can cause problems if 
there are some temperature dependent variables, beside 
those of the KLM, which can influence the magnetiza- 
tion. In our case this can be the JT splitting at special 



0(x) 



-E 



9e(k) 



5(x-e(k)) 



(CI) 



we now get 



+ 00 



4cosh^(i^) ' 



(C2) 

with the QDOS Paa{E, x) = -j^ImG^^iE) and the resis- 
tivity PciiT) — (T^^(T). Results with the same structure 
can be found in other worbiSiiSiSS. 
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